
Turbomachinery Laboratory 
Mechanical Engineering Department 


Texas A&M University 

College Station,Texas 77843-3123 






ENTRANCE AND EXIT REGION FRICTION FACTOR MODELS 
FOR ANNULAR SEAL ANALYSIS 


A Dissertation 
by 

DAVID ALAN ELROD 


Submitted to the Graduate College of 
Texas A&M University 

in partial fulfillment of the requirements for the degree of 
DOCTOR OF PHILOSOPHY 


August 1988 


Major subject: Mechanical Engineering 


ABSTRACT 


Entrance and Exit Region Friction Factor Models 
For Annular Seal Analysis. (August 1988) 

David Alan Elrod, B.S., Louisiana State University; 

M.S., Texas A&M University 

Co-Chairs of Advisory Committee: Dr. Dara Childs 

Dr. Clayton Nelson 

The Mach number definition and boundary conditions in Nelson’s nominally* 
centered, annular gas seal analysis are revised. A method is described for deter- 
mining the wall shear stress characteristics of an annular gas seal experimentally. 
Two friction factor models are developed for annular seal analysis: one model is 
based on fiat-plate flow theory; the other uses empirical entrance and exit region 
friction factors. The friction factor predictions of the models are compared to 
experimental results. Each friction model is used in an annular gas seal analysis. 
The seal characteristics predicted by the two seal analyses are compared to exper- 
imental results and to the predictions of Nelson’s analysis. The comparisons are 
for smooth-rotor seals with smooth- and honeycomb-stators. The comparisons 
show that the analysis which uses empirical entrance and exit region shear stress 
models predicts the static and stability characteristics of annular gas seals better 
than the other analyses. The analyses predict direct stiffness poorly. 
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NOMENCLATURE 


Direct and cross-coupled damping coefficients (FT /L) 
Diameter ( L ) 

Effective seal diameter, introduced in equation (22) ( L ) 
Honeycomb cell size, illustrated in figure 5 ( L ) 

Surface roughness (L) 

Seal reaction-force magnitude ( F) 

Fanning friction factor 
Radial clearance (L) 

Honeycomb cell depth, illustrated in figure 5 ( L ) 

Direct and cross-coupled stiffness coefficients (F/L) 
Seal length ( L ) 

Entrance length (L) 

Friction coefficients 

Seal inlet-to-exit pressure ratio (dimensionless) 

Fluid pressure (F/L?) 

Seal radius (L) 

Nominal duct-flow Reynolds number (l.ZZpUH / p) 
Flow-distance Reynolds number (pUx/p) 

Fluid bulk flow velocity ( L/T ) 

Dimensionless bulk flow velocity, U j Ru> 

Rotor to stator relative displacement components (L) 
Flow distance ( L ) 

Axial location within the seal (L) 


XI 


7 Ratio of specific heats for air 

k Entrance loss coefficient 

p Fluid viscosity (FT/L 2 ) 

p Fluid density (M/Z 3 ) 

t w Wall shear stress ( F/L 2 ) 

V Empirical friction multiplier, introduced in Chapter VI 

(dimensionless) 

u> Shaft angular velocity (1/T) 

Subscripts 

a Reservoir value 

b Sump value 

d Developing flow region 

r Rotor 

s Stator 

x Exit region 

z Axial 

6 Circumferential 

0, 1 Zeroth- and first-order values 
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CHAPTER I 
INTRODUCTION 

In pumps, compressors, turbines, and other turbomachines, annular seals are 
provided to limit the leakage between regions at different pressures. However, the 
forces within a seal may also influence rotordynamic response. In rotordynamics, 
a gas seal is modeled as a set of springs and dampers between the rotating shaft, 
or rotor, and the non-rotating housing, or stator. The model used to define the 
reaction-force/motion relationship for a nominally centered gas seal is 

;]{?}• <« 

In this equation, X and Y are the displacement components of the rotor relative 
to the stator, and Fx and Fy are the components of the reaction force acting on 
the rotor. The direct stiffness K , cross-coupled stiffness k, direct damping C, and 
cross-coupled damping c are referred to collectively as rotordynamic coefficients. 
The cross-coupled terms account for the fact that motion in one direction causes 
a force in an orthogonal direction. 

Figure 1 shows the forces on a rotor in a forward circular orbit of amplitude A 
and frequency w. Positive direct stiffness and cross-coupled damping act to center 
the rotor. Positive direct damping acts opposite to the velocity direction, opposing 
the whirling motion. However, positive cross-coupled stiffness acts to support the 
whirling motion-a destabilizing effect. Cross-coupled stiffness depends on the 
magnitude and direction (with respect to rotor rotation) of the circumferential 
velocity of the fluid within the seal. 


Journal Model: ASME Journal of Tribology 
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Since Lomakin (1958) first demonstrated that the stiffness of an annular, 
liquid seal can cause an increase in the critical speed of a rotor, a number of 
analyses have been developed to predict the rotordynamic coefficients of incom- 
pressible and compressible annular seals. For incompressible seals, Black and 
Jenssen (1969, 1970) use a bulk-flow analysis in which the circumferential bulk- 
flow velocity is assumed to be one half of the rotor surface speed and the wall 
shear stresses are defined by average circumferential and axial Reynolds numbers. 
Their “short-seal” solution neglects the effect of pressure perturbations on cir- 
cumferential flow. In another analysis, Black and Jenssen (1971) define wall shear 
stresses as a function of the local Reynolds numbers. Allaire et al. (1976) use 
Black’s model for small rotor motion about a centered position to numerically 
calculate dynamic coefficients for motion about an eccentric position. Black et al. 
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(1981) account for the influence of inlet swirl, demonstrating that the circumfer- 
ential bulk-flow velocity of a fluid element approaches one half of the rotor surface 
speed as it proceeds axially through the seal. Von Pragenau (1982) proposed a 
“damper seal” configuration, using a smooth rotor with a rough-surface stator, to 
reduce the fluid circumferential velocity and destabilizing cross-coupled stiffness of 
the seal. Childs’ (1983, 1984) incompressible seal analyses provide “finite-length” 
solutions for constant-clearance and tapered seals, in which the effect of pres- 
sure perturbations on the circ umf erential bulk-flow velocity are included. Childs 
and Kim (1985) have developed a finite-length solution for incompressible seals 
with different surface roughness conditions on the rotor and stator, the damper- 
seal configuration proposed by von Pragenau. Dietzen and Nordmann (1986) 
use finite-difference methods and a three-dimensional fluid flow model to calcu- 
late seal rotordynamic coefficients. Their analysis has been used for grooved and 
smooth-surface seals. 

Compressible flow in annular seals has been analyzed by Fleming (1978, 1980), 
Nelson (1984, 1985), and Nordmann et a). (1987). Fleming presents a short- 
seal solution for the leakage, direct stiffness, and direct damping of straight and 
tapered, smooth, annular gas seals. However, since he assumes one-dimensional, 
axial flow, Fleming cannot predict the cross-coupled coefficients. In his finite 
length analyses, Nelson accounts for the influence of inlet swirl and provides for 
shear and pressure effects on the circumferential bulk-flow velocity in constant- 
clearance or convergent-tapered gas seals with different rotor and stator surface 
roughnesses. In Chapter D, Nelson’s analysis is reviewed, and his Mach number 
definition and boundary conditions are revised. The revised analysis is referred 
to as the fully-developed-flow analysis in the present study. Nordmann et al. 
(1987) extend the finite-difference, incompressible seal analysis of Dietzen and 
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Nordmann (1986) to calculate rotordynamic coefficients for annular gas seals. 
For a honeycomb-stator seal, the stator-wall boundary conditions assumed by 
Nordmann et al. would not apply. 

For a given set of seal operating conditions, each of the above analyses predicts 
a pressure profile within the seal, which is integrated to determine forces and 
rotordynamic coefficients. Each analysis uses an abrupt loss at the seal entrance 
to model the high friction in the “developing flow region” of the seal. Actually, 
the developing flow region is often a substantial portion of the seal length. In 
Chapter ID, an improved model for friction in the entrance region of a seal is 
developed from flat plate flow theory. The flat-plate-flow friction model is used in 
a gas seal analysis (called the flat-plate-flow analysis), described in Chapter IU. In 
Chapter V, the friction factor predictions of the flat-plate-flow model are compared 
to experimental friction factors obtained using the test apparatus described in 
Chapter IV. In Chapter VI, experimental friction factor data are used to develop 
an empirical friction factor model for an annular gas seal analysis (the empirical- 
friction-factor analysis). In Chapters VII and VIII, static and dynamic results 
from experimental tests of four annular seals are compared to the predictions of 
the fully-developed-flow, flat-plate-flow, and empirical-friction-factor analyses. 



5 


CHAPTER n 

NELSON’S ANALYSIS: 

A FULLY-DEVELOPED-FLOW MODEL 

Nelson’s analyses (1984,1985) predict both static and dynamic characteristics 
for nominally-centered, annular gas seals. The static results include fluid leakage 
through the seal, pressure gradient along the seal axis, and the fluid axial and 
circumferential velocities through the seal. Dynamic data provided by the analyses 
consist of the rotordynamic coefficients (direct and cross-coupled stiffness and 
damping terms) in equation (l). 

In his analyses, Nelson uses a modified Hire’ (1973) turbulent bulk-flow fluid 
model to develop governing axial and circumferential momentum equations, and 
his model is completed by the continuity and energy equations. The equations 
define the relationship between the clearance, pressure, density, axial velocity, 
and circ umf erential velocity ( H , p, p, U 2 , and Ug) as functions of the spatial 
variables 6 and z and time t. Hire’ model uses the Blasius shear stress formula 
for fully-developed turbulent pipe flow, in which the friction factor / is 


/ = noR , 


mo 

e 


( 2 ) 


R e is the Reynolds number relative to the surface upon which the shear stress is 
acting, and the constants no and mo are generally determined empirically from 
static pressure flow experiments. In a bulk-flow analysis for incompressible seals, 
Nelson and Nguyen (1987) use Moody’s pipe friction formula (Massey, 1979) in- 
stead of equation (2). Previously used in seal analysis by von Pragenau (1982), 
this formula approximates duct friction as presented in the Moody chart: 


( z 10 6 \ 1/3 
= 0.001375 1 + 20000— + — 1 


) 


( 3 ) 
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where e is the surface roughness (“bump” height) of the duct, D is the effective 
diameter of the duct, and R e is the Reynolds number relative to the surface upon 
which the shear stress is acting. 

Ass umi ng small motion of the rotor about a centered position, Nelson uses a 
perturbation analysis similar to that employed by Childs (1983) to develop zeroth 
and first-order perturbation equations. The governing equations are expanded in 
the perturbation variables 


H = Ho + cHi p = po + epi p = po + tpi 

U x = Ugo + tUsi Ue = Ueo + tU$i. 


( 4 ) 


The zeroth-order solution represents a zero-eccentricity flow condition, with rotor 
rotation but without precession. This solution is developed iteratively and yields 
the mass-leakage flow rate, and the axial distributions of pressure, density, axial 
velocity, and circumferential velocity. The iterative solution scheme uses initial 
guesses for the zeroth-order seal entrance Mach number and entrance pressure-loss 
coefficient. The entrance-loss relationship is defined by 


Po(0) = 



(7-1) (« + 1) Mp (0) ' 
2 


t/U-'j) 


( 5 ) 


where p 0 (0) is the ratio of the seal entrance pressure to the reservoir stagnation 
pressure and Afo(0) is the entrance Mach number. The entrance Mach number 
is iteratively adusted, and the loss coefficient k is recalculated according to the 
relationship 


k = 


5.3 


log 10 Re 


- 1.0 . 


( 6 ) 


k is used to model the high friction predicted by Deissler (1953) in the developing 
flow region of a duct. The iterative solution procedure for Mo(0) and k continues 
until either: 
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1) the Mach number at the exit reaches unity and the exit pressure is 
greater than the sump pressure (choked flow) , or 

2) the exit pressure equals the sump pressure and the exit Mach number 
is less than unity (unchoked flow). 

The pressure, density, and velocity distributions and their axial derivatives, 
which are determined in the zeroth-order solution, are used in defining coefficients 
of the first-order perturbation equations. These equations define the pressure, 
density, and axial and circumferential velocity perturbations (pi, pi, U a i, and 
U$ i) due to rotor motion. The four physical boundary conditions required for 
the solution of these equations depend on the perturbation conditions that are 
specified at the seal entrance and exit. Numerical integration of the governing or- 
dinary differential equations in the perturbation variables provides the first-order 
solution. The first-order pressure solution is then integrated axially and circum- 
ferentially for the rotordynamic coefficients in equation (l). The seal boundary 
conditions used by Nelson (1984, 1985) have been changed for the fully-developed- 
flow, flat-plate-flow, and empirical-friction-factor analyses of the present study. A 
discussion of the changes follows. 

In Nelson’s analyses (1984, 1985), the zeroth-order circumferential velocity at 
the seal entrance (i.e. at z = 0) is equal to the reservoir circumferential velocity, 
and the entrance circumferential velocity perturbation is zero; i.e., 

U 0 O {O) = U O a U 0l { 0)=0. (7) 

These boundary conditions have not been changed in the present analyses. 

The Mach number definition in Nelson’s analyses is 

= ^ . 

IV 


( 8 ) 
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For choked flow, M = 1 at the seal exit. When Mach number is defined by 
equation (8), the fluid velocity exceeds the velocity of sound when M = 1 unless 
U$ = 0. In general, the circumferential component of the fluid velocity is not zero 
at the seal exit. In the present study, the Mach number in the seal is defined by 


M 2 (ul±uj)p 

IP 


(9) 


The fluid velocity equals the velocity of sound when M = 1 in equation (9). The 
first-order perturbation of the Mach number in equation (8) is 


= ■ (10) 
For choked flow, M\ is zero at the seal exit (z = L). In Nelson’s analyses, this 


yields 


Pi(L) , V.i(L) Pi(L) 


+ 


= 0 


( 11 ) 


2po(L) U z o(L) 2po(L ) 

Using the Mach number definition of equation (9), equation (ll) is replaced by 


Pi{L) U z o(L)^iW + U eo (L)U gl (L) Pl {L) 

2/> 0 (X) %(I) + Uj 0 {L) 2po (L) 


For unchoked flow, the first-order perturbation in the exit pressure is zero: 


Pi(L) = 0 . 


(13) 


An expansion of equation (5) in the perturbation pressure and perturbation 
Mach number (equation (10)) yields a first-order pressure loss equation which 
must be satisfied at the seal entrance in Nelson’s first-order solution: 


(n) 4. r Po(°) ( ElB 4. EnB \ 

PUUJ+ 2p 0 (0)-r VPo(0) + U a0 (0)J 


= 0 , 


(14) 


where 


r = 


7(/c-H)po(0)Mg(0) 


l + b-l) (»c + l) 


A^(0) * 


(15) 


2 
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In the present analyses, equation (14) is replaced by 

„ rm x gg(g> /' Mg) x 2C W°)^,i(o) \ _ n 

PlU 2Po(0) - T U(O) %(0) + %(0) ; ~ 

which is the first-order expansion of equation (5) using equation (9), with t/ei(0) = 

0. 



The final first-order boundary condition in Nelson’s analyses is derived from 
an expression for the density at the seal entrance: 


Po(0) = Po(0) + 



(17) 


Expanding equation (17) in terms of the perturbation pressure, density, and Mach 
number (equation (10)), Nelson obtains 


Px(0) + 


r *Po(0) 
2p 0 (0) + rs 


( 2U, i(0) 

V ^o(O) 



(18) 


where 


a = 1 + (Tf - 1) 


(Afo(0) («+!) — 2] 


(19) 


27 (/c + 1) 

and r is defined by equation (15). In the present analyses, use of the Mach number 
definition in equation (9) alters this boundary condition as follows: 


?.(«) ~ lw ?l(0) + ~ 1>?0 7 ( ° >t '' ofo) t, .i(°) = 0 • 


( 20 ) 


In the present study, the boundary conditions in equations (7), (12), (16), and 
(20) are used in the fully-developed-flow, flat-plate-flow, and empirical-friction - 
factor analyses, with k = 0 for the flat-plate-flow and empirical-friction-factor 
analyses. 
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CHAPTER HI 

A FLAT-PLATE-FLOW 
FRICTION FACTOR MODEL 

The “developing flow” or “entrance” region of a seal is often a substantial 
portion of the seal length. However, in seal analyses, an entrance-loss coefficient 
(e.g. k in equations (5) and (6)) is generally used to model the high losses in 
the developing flow region as abrupt losses at the seal inlet, and fully developed 
flow is assumed throughout the seal. In this chapter, a friction factor model is 
developed for the entrance region of an annular seal as the following questions are 
answered: 

1) In duct flow, what are the developing flow and fully developed flow 
regions? 

2) How does fluid friction in the developing flow region differ from fluid 
friction in the fully developed flow region? 

3) In the flat-plate-flow seal analysis, how does the friction model in 
the developing flow region differ from that of the fully-developed-flow 
analysis? 

As fluid enters a duct, a viscous boundary layer is formed along the duct wall 
(figure 2). At the wall, the fluid velocity is zero. Viscous effects cause a reduction 
of the velocity of adjacent fluid particles. At the center of the duct, a “free stream” 
region exists where the wall shear has not yet caused a reduction in fluid velocity. 
The region between the free-stream core and the duct wall is the boundary layer. 
The boundary layer becomes thicker as the fluid flows through the duct, and the 
free stream accelerates, satisfying the conservation of mass requirement. In figure 
2, a plot of the velocity u versus the radial location r at a given z is called the 
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velocity profile at z. If the duct is long enough, the boundary layer will extend 
to the center of the duct, and at some finite distance L e (the “entrance length”) 
from the duct entrance, the velocity profile will cease to change with increasing z. 
When the velocity profile ceases to change, the flow is said to be fully developed. 
The entrance or developing flow region extends from the duct inlet to the fully 
developed flow region. 



Figure 2. Developing flow in a duct. 


Fluid friction models for the fully developed flow region are numerous (see 
Serghides, 1984). The accuracy of a model is generally evaluated by comparing 
the friction factor (/) predicted by the model to friction data plotted in the Moody 
diagram. The relationship between friction factor / and wall shear stress r w is: 

r Tw 

J ~ 0.5pU 2 ' 

In their seal analyses, Nelson (1984) and Nelson and Nguyen (1987) use two fully 
developed flow friction models: Moody’s equation, which approximates the Moody 
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diagram well for Reynolds numbers between 4000 and 10 7 and relative roughnesses 
(e/D) up to 0.01, and Blasius’ equation. For a given pipe Reynolds number, 
Blasius’ model can be used to approximate both the Moody diagram friction factor 
/ and the effect of small changes in Reynolds number on /. In general, however, 
the friction coefficient values no and mo which must be used to approximate 
the Moody diagram depend on Reynolds number and relative roughness. In the 
friction factor model developed in this Chapter, Moody’s equation is used for the 
fully-developed-flow friction factor when the relative roughness is less than 0.01 
("smooth’’ surfaces) and Blasius’ model is used for honeycomb surfaces. 

Deissler (1953), Shapiro and Smith (1948), and Keenan and Neumann (1946) 
are among those who have investigated fluid friction in duct entrance regions. In 
an analysis of turbulent flow in the entrance regions of smooth passages, Deissler 
(1953) extends a previous analysis for fully developed, turbulent flow. Figures 
3 and 4 (from Deissler, 1953) show the friction factors Deissler predicts for flow 
through a tube and for flow between parallel plates, respectively. In each plot, 
X/D is the flow distance in hydraulic diameters. Deissler predicts that the friction 
factor reaches its fully developed flow value in about 10 diameters. In his annular 
seal analysis, Nelson (1984) uses a fully developed flow friction factor throughout 
the seal, modeling the excess losses predicted by Deissler in the entrance region 
as abrupt losses at the seal inlet. 

In an experimental study in the 1940’s, Shapiro and Smith (1948) determined 
friction factors in the entrance region of smooth, round tubes with bellmouth 
entries for uniform inlet conditions. Shapiro and Smith used water and air at low 
Mach numbers as test fluids, at Reynolds numbers (based on tube diameter) from 
39,000 to 590,000. The diameters of the test pipes ranged from 9.5 to 102 mm. 
Their study shows that, for developing duct flow, the Reynolds number based on 
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distance from the inlet is more significant in predicting friction factors than the 
pipe (diameter-based) Reynolds number. The friction factor is within 5% of the 
fully developed flow value “beginning at about 40 to 60 tube diameters from the 
tube inlet.” Furthermore, Shapiro and Smith predict that the friction factor is 
“almost infinite” at the tube inlet. Boundary layer characteristics in the entrance 
region of a tube are similar to boundary layer characteristics in flat plate flow 
studies. 

In another experimental study, Keenan and Neumann (1046) report friction 
factors for subsonic air flow at Reynolds numbers from 100,000 to 450,000. The 
test pipe was a piece of standard drawn-brass tubing with an inside diameter of 
9.5 mm. The study shows that the friction factor decreases from a high value near 
the tube inlet to the fully developed flow value at a flow distance of about 40 pipe 
diameters. 

For friction-factor calculations, the effective diameter of an annular seal is 
four-thirds of the seal radial clearance (White, 1979); i.e., 

D' = l B . (22) 

Based on the data described in the preceding paragraph, fully developed friction 
factor values may not be attained for a substantial portion of the seal length. 
For example, forty diameters is over 60% of the length of the turbine interstage 
seal of the High-Pressure Oxidizer Turbopump of the Space Shuttle Main Engine 
(described by Nelson, 1984). Previous seal analyses (and the fully-developed- 
flow analysis of the present study) use fully-developed-flow friction models, which 
are a function of duct diameter, throughout the seal. The only exception is the 
analysis of Mullins and Barrett (1983), which uses an entrance region friction 
factor model. However, the friction factor is assumed to increase from zero at 
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the seal inlet to the fully-developed-flow friction value at the end of the entrance 
region, which disagrees with the results described in the previous paragraphs and 
shown in figures 3 and 4. Furthermore, the entrance region friction factor used 
by Mullins and Barrett depends on the effective diameter of the seal and on the 
axial component of the flow distance. In the following paragraphs, a new friction 
factor model for the entrance region of a duct is developed. In this model, friction 
is not a function of duct diameter; it is a function of the total flow distance within 
the seal (axial and circumferential components). 

Fluid friction in the entrance region of a duct has been compared to fluid 
friction along a flat plate. For flat plate flow (zero pressure gradient by definition) , 
/ depends on distance from the leading edge of the plate. A number of equations 
have been used to model / for turbulent flow past a flat plate, and a common 
functional form is 

/ = aRj , (23) 

where x is the distance from the leading edge and is the Reynolds number 
defined by: 

n 

If a pressure gradient is present (i.e. dp/dx ^ 0), / is determined by numerical 
analysis. According to White (1974), there are few analyses for compressible, 
turbulent flow with a non-zero pressure gradient over a flat plate. Available solu- 
tion methods require knowledge of the free-stream velocity gradient and empirical 
boundary layer relationships. The use of one of these plate-flow numerical solution 
methods (an external flow analysis) for entrance flow in a duct (an internal flow 
problem) is questionable since boundary layer interaction in a duct complicates 
the flow. 
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In the flat-plate-flow analysis of the present study, Nelson’s (1984, 1985) 
governing equations are used, and the solution method is the same as for the 
fully-developed-flow model. In the flat-plate-flow analysis, the friction factor in 
the developing flow region of an annular seal is modeled by equations (23) and (24) , 
where i is a flow distance relative to the rotor or stator surface. For calculating 
the wall shear stress at the seal stator in the developing flow region, the friction 
factor for the zeroth-order solution is 

U = «. {~) > ” Wh + tf, 2 o) IW,) • (25) 

For the shear stress at the rotor, the friction factor is 

= . (26) 

The equations for flow distance relative to the stator and rotor ( x t and x r , re- 
spectively) at an axial location Z are 



The coefficient a and the exponent b in equation (23) are calculated by setting the 
developing-flow friction factor /<* equal to the fully-developed-flow friction factor 
// at a flow distance of forty diameters; i.e., 

/d(*=40D,) = o = fl(40/Ee) b = // . (29) 

When Moody’s equation (equation (3)) is used for fj, a and b values are calculated 
by solving the following equations: 

/, = o(40He) fc = 0.001375 fl + (20000^-jlg- + ' 1 (30) 
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and 


|§-=#(40.R«)‘ = 1375 

dlt. R , ' cl 


3 Rl 


( 


t 10 6 \ 2/3 

20000 —— + — 

1.33 H 0 J2e ) 


(31) 


For any iteration in the zeroth-order solution (described in Chapter II), only a 
and b are unknown in equations (30) and (31). Equations (30) and (31) are used 
to calculate a set of coefficients for both the rotor and the stator (i.e. a, , b r and 
a t , b t ). 

The expansion of equations (25) and (26) in terms of the perturbation vari- 
ables in equation (4) yields 


f.i — f»obt 

frl = /rO&r 


Pi UeoUoi + U*V,1 

.Po XiO UZo + U?o . 

Pi Zri_ (Ueo — Rv)Uei + U^JJzi 

.Po *ro (Ue 0 - Ru) 2 + Ul 0 


(32) 

(33) 


which are required for the first-order solution in the flat-plate-flow analysis. In 
equations (32) and (33), 


dXm dx s 

= euL u “ + eu; o u " 

(34) 

OXr OXr 

*'' = au m v ’' + au, a v ‘' ■ 

(35) 


Using Leibnitz’s rule and equation (27), 


dx t 

dU e o 





(36) 


Similar expressions can be derived for the other partial derivatives in equations 
(34) and (35). For the flat-plate-flow seal analysis, the following equations are nu- 
merically integrated in the zeroth-order solution for use in calculating the zeroth- 
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and first-order wall shear stresses in the seal: 


dx t 

17 

dx r 

dz 

a{dx,/dUeo) 

dz 

aidzjdiu) 

dz 

d{dxr/dUe o) 

dz 

a&xriav*) 

dz 



(37) 

(38) 

(39) 

(40) 

(41) 

(42) 


For a seal with a honeycomb stator, flow at the stator surface differs from the 


developing flow described above. In figure 5, the honeycomb cell size is d c and the 
cell depth is h c . The cell size and depth in a honeycomb-stator seal are usually 
from one to several times the seal clearance. The fluid velocity along the stator 
wall is not necessarily zero. Use of the three-dimensional seal analysis of Nord- 
mann et al. (1987), which assumes “no-slip” boundary conditions at the stator, is 
questionable for a honeycomb seal. In the flat-plate-flow analysis, the developing- 
flow friction factor defined above is only used for shear stresses at the rotor for 
a honey comb-stator/smooth-rotor seal. Throughout the seal, honeycomb-stator 
wall shear stresses are calculated using the Blasius friction factor: 


/ = noR™ 


( 2 ) 


Chapter V includes comparisons of the friction factor predictions of the above 


model to experimental friction factors. 
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Figure 5. Honeycomb stator insert detail. 
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CHAPTER IV 
TEST APPARATUS 

A complete description of the test apparatus used in this study is provided by 
Childs et al. (1986 and 1988). The rotor shaft is suspended, pendulum fashion, 
from an upper, rigidly-mounted, pivot shaft, as illustrated in figure 6. This ar- 
rangement allows for a horizontal (harmonic) motion of the rotor. A cam within 
the pivot shaft allows vertical (static) positioning of the rotor. The rotor is ex- 
cited, horizontally, by a hydraulic-shaker head which acts on the rotor-shaft hous- 
ing. The design of the test rig, illustrated in figure 7, permits the installation 
of various rotor/stator combinations. The stator is supported in the test-section 
housing by three piezo-electric, quartz, load cells in a trihedral configuration. 



Figure 6. Test apparatus. 
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Hydraulically Filled 



Figure 7. Test-section cross section. 


For the present study, the test apparatus has been used to determine fric- 
tion factors for the six seals listed in Table 1, and dynamic characteristics for 
seals SI, HCl, HC2, and EC3 of Table 1. In Table 1, seals SI, S2, and S3 are 
smooth-stator /smooth-rotor seals, and seals HCl, HC2, and HC3 are honeycomb- 
stator /smooth-rotor seals. In the honeycomb stator insert detail shown in figure 
5, the honeycomb cell size is d c and the cell depth is h c . In Table 1, the cell sizes 
are 1.57 mm and 0.51 mm, and the cell depths are 1.91 mm and 1.47 mm. In all 
six seals, the diameter and measured surface roughness of the rotor were 151.4 
mm and 0.188 /zm, respectively. 

The operator of the dynamic-seal-apparatus can control the following four 
independent variables: inlet pressure, rotor speed, shake frequency, and inlet 
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circumferential velocity. For dynamic tests, the actual test points for three of 
these variables are shown in Table 2. The sump pressure Pb was about 1.01 bars 
(atmospheric pressure) . For all dynamic test results in the present study, the rotor 
was shaken about a centered position within the stator at a shake frequency of 
74.6 Hz, and a shake amplitude of about 80 pm. 


Table 1. Seal Dimensions 


Seal 

Length 

Clearance, E 

Stator Surface 

SI 

5.08 cm 

0.41 mm 

e = 0.188/zm 

S2 

5.08 cm 

0.30 mm 

e = 0.188/xm 

S3 

5.08 cm 

0.20 mm 

e = 0.188/xm 

HC1 

5.08 cm 

0.41 mm 

1.57x1.91 mm cells 

HC2 

5.08 cm 

0.41 mm 

1.57x1.47 mm cells 

HC3 

5.08 cm 

0.41 mm 

0.51x1.47 mm cells 


Table 2. Test Variables 


Inlet Pressure 

Rotor Speeds 

Inlet Circumferential Velocities 

1- 3.08 bars 

2- 4.46 bars 

3- 5.84 bars 

4- 7.22 bars 

5- 8.26 bars 

1- 3000 cpm 

2- 6000 cpm 

3- 9500 cpm 

4- 13000 cpm 

5- 16000 cpm 

-2-High velocity against rotation 
-1-Low velocity against rotation 

0- Zero circumferential velocity 

1- Low velocity with rotation 

2- High velocity with rotation 


The inlet circumferential velocity is controlled using four sets of inlet guide 
vanes. For the two “high” inlet circumferential velocity cases in Table 2, the 
guide- vane depth (“A” in figure 8) is 35% of the guide- vane depth for the “low” 
inlet circ umf erential velocity cases. The zero inlet circumferential velocity case 
is obtained without guide vanes. For a set of swirl vanes at a constant running 
speed, inlet circumferential velocity remains almost constant for the test pressure 
ratios in the table above. The velocity decreases slightly with rotor speed, mainly 






Figure 8. Inlet-guide-vane detail. 


24 


because the rotor grows with increasing speed and reduces the leakage. The ratio 
of inlet circumferential velocity to rotor surface velocity, u$ a , ranged from about 
—5 to about +5. Average experimental inlet circumferential velocities for seals 
Si, HCl, HC2, and HC3 are given in Table 3. 


Table 3. Inlet Circumferential Velocities 


Swirl 

Smooth Seal 

HCl 

HC2 

HC3 

-2 

-100 m/sec 

-58 m/sec 

-68 m/sec 

-53 m/sec 

-1 

- 44 m/sec 

-26 m/sec 

-31 m/sec 

-24 m/sec 

0 

0 m/sec 

0 m/sec 

0 m/sec 

0 m/sec 

1 

45 m/sec 

27 m/sec 

31 m/sec 

24 m/sec 

2 

118 m/sec 

68 m/sec 

80 m/sec 

62 m/sec 




CHAPTER V 


EXPERIMENTAL FRICTION FACTORS 

The test apparatus described in Chapter IV has been used to obtain friction 
factor data for the six seals listed in Table 1 of Chapter IV. In the present Chapter, 
the experimental method is described, and friction factor test results are compared 
to the predictions of the friction model used in the flat-plate-flow seal analysis. 
Test Method 

One-dimensional, steady, adiabatic flow of a perfect gas with constant specific 
heats through a constant-area duct with no external work is called Fanno line flow. 
The momentum equation for a control volume of length dz is 


-Adp - t w A w = pAUdU , 


(43) 


where A is the cross-sectional area of the duct, r w is the shear stress at the duct 
walls, and A w is the surface area over which t w acts. For a seal with the rotor 
centered, 


A = - (Dj ~ D ,) , A w = n(D e + D r )dz , and D, - D r = 2 H 0 . (44) 

4 


Substitution and simplification yields 

— dp — ~~~dz — pU dU . (45) 

no 

The relationship between friction factor / and r w is expressed in equation (21): 


/ = 


0.5pU 2 ' 


( 21 ) 


Using equations (21) and (45), the perfect gas law, the definition of Mach number, 
and the conservation of mass requirement, one can derive (see John, 1984) that 
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the effective friction factor, / e , for a seal is: 

= — 1-M 2 1 dM 

e “ M + M 2 )7M 2 J dz 

To determine f e for a seal using the test apparatus described in Chapter IV, 
the stator is mounted in the apparatus, a smooth rotor is centered in the stator, 
and flow tests are conducted. For each test, the rotor remains stationary (no 
rotation), the inlet air has no circumferential velocity component, and the seal 
inlet pressure is held constant. The stagnation temperature of the air, flow rate 
through the seal, and pressures along the seal length are measured. The Mach 
number is calculated at 3.2 cm (0.125 inch) intervals along the seal (i.e. at pressure 
tap locations) . A least-squares cubic curve fit is used to express the Mach number 
as a function of axial location. Equation (46) is used to calculate f e in the seal. 
Values of f e have been calculated for three smooth-stator/smooth-rotor seals (SI, 
S2, and S3 in Table 1), and for three honeycomb-stator/smooth-rotor seals (HCl, 
HC2, and HC3 in Table 1). 

For a seal with identical rotor and stator surfaces (i.e. a smooth-stator /smooth- 
rotor seal), 

/. = Jr = /. . (47) 

For a seal with different surface roughnesses (e.g. a honeycomb-stator /smooth- 
rotor seal), / e , /,, and f r are related as follows. In equation (43), 

t w Axi/ — t w (A* -f- Af) = 7|i4j "h t t A t , (48) 



and 

A = («) 

for a seal. Substitution from equation (21) and (49) into equation (48) yields 
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Rearranging, 

/. = (l + ^)/.-^/r • (51) 

To determine the friction coefficients no, and mo, for a honeycomb stator, 
stationary-rotor tests are run, and f, values are calculated for nine inlet pressures 
using equations (46) and (51). In equation (51), f r is calculated using smooth- 
stator /smooth-rotor seal test results. When friction coefficients are calculated, 
entrance-region and end effects are avoided by using test data for axial locations 
from forty effective seal diameters to about 60% of the seal length. For each stator, 
average f, results and average experimental Reynolds numbers in this region are 
used in a least-squares curve fit to determine mo, and no, (see equation (2)). The 
friction coefficients mo and no for the six seals listed in Table 1 are given in Table 
4. 


Table 4. Seal Friction Coefficients 


Surface 

no 

mo 

Smooth 

0.154 

-0.299 

HCl stator 

0.303 

-0.117 

HC2 stator 

0.155 

-0.098 

HC3 stator 

0.424 

-0.134 


Experiment /Theory Comparisons 

For the six seals of Table 1, typical plots of f e versus dimensionless axial 
location ( z/L ) are shown in figure 9. For seal S2, the test case shown is for 
unchoked flow. For the other plots in figure 9, flow is choked. In each plot, 
experimental f e results are designated by the letter “x" . The error bands represent 
the uncertainty determined using the method described by Holman (1978). The 




FRICTION FRCTOR 
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solid lines in the plots represent the friction factors predicted by the flat-plate-flow 
friction model and equations (47) and (50). For seal Si of Table 1, 40 D e « 0.4L. 
For seal S2, 40 D e « 0.3Z,. In figure 9, the experimental friction factors for seals 
SI and S2 are high and decreasing for about forty diameters. For the honeycomb- 
stator seals, f e decreases in the first one- third of the seal, and then increases. For 
all six seals, the experimental friction factor is generally higher than predicted 
in the entrance region. In the middle 20% to 30% of seals Si and S2, f e is 
higher than predicted by Moody’s equation. In the last 40% of all six seals of 
Table 1, f e decreases experimentally. No existing friction factor model, including 
the flat-plate-flow friction model, predicts this decrease in /. In Chapter VI, a 
friction factor model is developed from experimental data obtained using the test 
apparatus described in Chapter IV. 
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CHAPTER VI 

AN EMPIRICAL FRICTION FACTOR MODEL 

Figure 9 in Chapter V shows that the friction factor in a seal decreases in the 
“exit" region of the seal. No existing friction factor model includes this feature. 
The high friction in the “entrance" region of the plots in figure 9 is underpredicted 
by the flat-plate-flow friction model. Furthermore, Moody’s equation underpre- 
dicts the friction factor between the entrance and exit regions of the smooth- 
stator /smooth-rotor seals. In this chapter, a friction factor model based on the 
data in figure 9 is developed for annular seal analysis (the “empirical-friction- 
factor" seal analysis). The friction factor predictions of the model are compared 
to experimental data from stationary- rotor tests of the six seals of Table 1. In the 
empirical-friction-factor seal analysis, the friction model in this chapter is used in 
Nelson’s (1984,1985) governing equations for a gas seal, and the solution method 
is the same as for the fully-developed-flow model. 

For calculating the wall shear stress at the seal stator, the zeroth-order friction 
factor is 

/• o = V<uoVxofB«o • (52) 

For the shear stress at the rotor, the friction factor is 

frO = VdroVxofBrO • (53) 

In equations (52) and (53), and ford are the zeroth-order Blasius friction 
factors for the stator and rotor, respectively; i.e., 

/b.o = no. [Via + (54) 

/BrO = no. {(U K - Ru)* + 


(55) 
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The terms ¥4, o, ¥4 r o, and ¥*o in equations (52) and (53) are defined in the 
following paragraphs. 

Exit Region friction 

No existing friction model predicts a decrease in friction factor near the seal 
exit, as shown in figure 9. Others (e.g., Keenan and Neumann, 1946) have rec- 
ognized that “end effects” are present in duct flow, but have chosen to avoid the 
region, experimentally. Figure 9 shows that the friction factor decreases to near 
zero a s z/L increases from 0.6 to 1.0. In the empirical-friction-factor seal analysis, 
this region (0.6 < z/L < 1.0) is defined as the “exit region” of the seal, and ¥*0 
is used in equations (52) and (53) to model the friction factor decrease, where 

¥*0 = -4.12 -I- 20.3 (z/L) - 25.1 (z/L) 2 + 9.03 [z/L) 5 . (56) 

when z/L > 0.6. When z/L < 0.6, ¥*o = 1.0. The constant terms in equation 
(56) were obtained using the following friction factor data: 

= 0.60) =0.00 
¥* 0 [z/L = 0.60) = 1.00 
¥ x0 (z/L = 0.80) =0.72 
V x0 {z/L = 1.02) =0.00 

Entrance Region Friction - Smooth Surfaces 

In equations (52) and (53), ¥4* and ¥a r ere used to model the high friction 
in the developing flow region, where 

¥* = 2.76 - 5.1 6(x t /L e ) + 5.04(x t /L e ) 2 - lM(xjL c f (57) 

and 

¥dr = 2.76 - 5.16 {x r /L e ) + 5.04(x r /L e ) 2 - 1.64(x r /X e ) 3 , (58) 
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when x/L e < 1.0, and 

L e = 40 D e = 53.33# . (59) 

For x/L e > 1.0, = 1.0. The constants in equations (57) and (58) were obtained 

using the following friction factor data: 

1.00) =0.00 
* d {x/L e = 0.14) = 2.13 
* d {x/L e = 0.47) = 1.28 
* d {x/L e = 1.00) = 1.00 

Figure 10 shows a comparison of typical experimental and predicted friction fac- 
tors for the smooth-stator seals of Table 1. The model used for the predictions in 
f - re 10 is defined by equations (52), (54), (56), and (57). 

Honeycomb Stator Friction Factors 

For all honeycomb-stator seals, equation (57), which is used in the calculation 
of smooth-stator friction factors, is replaced in the present empirical friction factor 
model by 

= 5.68 - 36.5 (z/L) + 83.5 (z/L) 2 - 58.4 (z/L) 3 (60) 

for z/L < 0.6. For z/L > 0.6, = ^^(z/L — 0.6). The constants in equation 

(60) were obtained using the following friction factor data: 

. 0-60) =0 
¥*(*/£ = 0.06) = 3.7 
* dt (z/L = 0.33) = 0.63 
*ds(z/L = 0.50) = 1.00 

Figure 11 shows a comparison of experimental and predicted f e for seals HCl, 
HC2, and HC3 of Table 1. The model used for the predictions in figure 11 is 
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Figure 10. Empirical vs. experimental f e for the smooth-stator seals of 
Table 1. Solid lines: empirical predictions. 
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defined by equations (50), (52) through (56), (58), and (60). The experimental 
results show that f e is high near the inlet of the honeycomb-stator seals, which 
disagrees with the friction factor model used in the flat-plate-flow analysis. 
First-Order Perturbations 

For the empirical friction factor model above, the first-order perturbations of 
equations (52) and (53) in terms of the variables defined in equation (4) are 


and 


f. 1 = /. 0 ( 


*J.l [ *.I /b 

ViU 0 ¥*0 f B 


h. i\ 

fBtO/ 


Ui = u ( 


*drl | [Brl\ 

9dr0 0 /Bi-0 / 


(61) 


(62) 


Since is independent of the perturbation variables, ^*1 = 0. Similarly, = 
0 for a honeycomb stator. For a smooth stator, 





Ei 

Ho 


for x i0 /L e < 1.0, where 


rl>, = -5.16(x, 0 /L e ) + 10.08(x, 0 !L e ) 2 - 4.92(x, 0 /I e ) 3 . 


and x t i is defined by equation (34). For x l0 > 1-0, = 0. 9^1 is expressed 

similarly. Finally, 


/b» i = fB»o^io a 


Pi H± UeoUjn + UzpUzi 

.po n o+^o 




and 


SBt\ — fBr0*nO T 


Pi (t/so — i?a>)t/gi -f UsqUz\ 

.po #o (tfso - Ru) 2 + Ul 0 
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CHAPTER VH 

STATIC PRESSURE AND LEAKAGE RESULTS 

The test apparatus described in Chapter IV has been used for stationary' 
rotor tests of six seals. Three smooth-stator/smooth-rotor seals (SI, S2, and S3 
in Table l) and three honey comb-stator /smooth- rotor seals (HCl, HC2, and HC3 
in Table 1) have been tested, with no prerotation of the inlet air. Stationary-rotor 
tests were conducted at nine inlet pressures: from 3.08 bars to 8.60 bars, in 0.69 
bar increments. Seals SI, HCl, HC2, and HC3 of Table 1 have also been tested 
dynamically; the dynamic test conditions are described in Chapter IV. 

In this chapter, experimental leakage and pressure gradient results are com- 
pared to the predictions of three gas seal analyses. In the figures in this chapter, 
curves I, 2, and 3 represent the predictions of the fully-developed-fiow, fiat-plate- 
Sow, and empirical-friction-factor analyses, respectively. Experimental results are 
designated by the letter u x”. The experimental results have been reported by 
Elrod and Childs (1988). 

Leakage 

All three analyses predict leakage well, especially for the honeycomb seals. 
Typical plots of predicted and experimental leakage for the test seals are shown in 
figure 12. The results shown are from dynamic tests of seals Si and HC2 of Table 
1. The uncertainty in the experimented results is about 1%. In the plot for seal 
Si, the circumferential velocity ratio is the ratio of the circumferential velocity of 
the air upstream of the seal to the rotor surface speed [Ue a /Ru). The empirical- 
friction-factor analysis predicts lower leakage values than the fully-developed-flow 
and flat-plate-flow analyses because the friction factor in the entrance region of 
the seal is much higher for the empirical-friction-factor analysis. 


LEAKAGE (kg/sec) LEAKRGE (kg/sec) 


SEAL 51 , P a *=B . 2 E bars 
ROTOR SPEED-13000 cpm 



5 ERL HC2 

SWIRL: 0, Pa=3.0B bars 



ROTOR SPEED (cpm) 


Figure 12. Leakage results for dynamic tests of seals Si and HC2 of 
Table 1. 



Pressure Gradient 

Experimental and predicted pressure gradients are shown in figure 13. The 
results shown are from dynamic tests of seals Si and HC2 of Table 1. Both test 
cases shown in the figure are choked flow cases. P/Pa is the dimensionless seal 
pressure, and z/L is the dimensionsless axial location. The uncertainty in the 
dimensionless pressures is less them 0.05. In each test case, there is little dif- 
ference between the pressure gradients predicted by the fully-developed-flow and 
flat-plate-flow analyses. The flat-plate-flow and empirical-friction-factor analyses 
predict a steeper slope for the pressure gradient at the seal inlet (z/ L< 0.1) than 
the fully-developed-flow analysis, because the difference between the friction fac- 
tors of the flat-plate-fiow and empirical-friction-factor analyses and the friction 
factor used in the fully-developed-flow analysis is greatest at the seal inlet. 

Generally, the empirical-friction-factor analysis predicts the pressure gradient 
better than the fully-developed-flow and flat-plate-flow analyses, especially in the 
exit region (z/L > 0.6) of the seals. The pressure drop in the last 10% to 25% 
of the seal length is overpredicted by the fully-developed-flow and flat-plate-flow 
analyses. For seal HC2, the pressure gradient predicted by the empirical-friction- 
factor analysis for z/L < 0.4 indicates a need for further investigation of the 
friction factor characteristics of honeycomb surfaces. 
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CHAPTER Vm 
DYNAMIC TEST RESULTS 

The air-seal apparatus described in Chapter IV has been used to test 
smooth-stator/smooth-rotor seals with constant and tapered stator diame- 
ters, honeycomb-stator seals with smooth and labyrinth rotors, labyrinth- 
stator /smooth-rotor seals, smooth-stator/labyrinth-rotor seals, and a seal with 
interlocking stator and rotor. In this chapter, rotordynamic coefficient results 
from tests of one smooth-stator/smooth-rotor seal (SI in Table 1, Chapter IV) 
and three honey comb-stator /smooth-rotor seals (HCl, HC2, and HC3 in Table 
1) are compared to the predictions of three gas seal analyses. The predicted de- 
pendencies of rotordynamic coefficients on rotor speed (w), inlet circumferential 
velocity ratio (u$ 0 = U$ a /Rw) t and seal inlet-to-exit pressure ratio ( Pr = p a /pi ) 
are compared to experimental results. 

As described in Chapter I and illustrated in figure 1, the stability of a seal 
in which the rotor is whirling at a frequency u depends on the cross-coupled 
stiffness and direct damping ( k and C, respectively) of the seal. Cross-coupled 
stiffness depends on the magnitude and direction (with respect to rotor rotation) 
oiU e - In the following paragraphs, comparisons of the experimental and predicted 
dependencies of k and C on u t u$ a , and Pr are emphasized. Comparisons for the 
direct stiffness, K, are also presented. 

For the rotor-speed dependence of the rotordynamic coefficients, results are 
shown for ug a = 0 (swirl case 0 of Table 2). For any other swirl case, the inlet 
circumferential velocity ratio u$ a decreases by about 80% when the rotor 6peed is 
increased from 3000 cpm to 16000 cpm. For the dependence of the rotordynamic 
coefficients on the inlet circumferential velocity ratio, results are shown for 0 < 
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u$ 0 < 1.2. For the dependence of the rotordynamic coefficients on Pr , results 
are shown for u$ 0 « 1.0, a typical value of u$ a in practice. In the Bgures in this 
chapter, curves 1 , 2, and 3 represent the predictions of the fully-developed-Bow, 
Bat-plate- Bow, and empirical-friction-factor analyses, respectively. Experimental 
results are designated by the letter “x The experimental results have been 
reported by Elrod and Childs (1988). 

Relative Uncertainty 

The uncertainty in the dynamic coefficients can be determined using the 
method described by Holman (1978). The uncertainty in the force, excitation 
frequency, and displacement measurements are 0.44 N (0.1 lb), 0.065 Hz, and 
0.0013 mm (0.05 mils), respectively. The maximum calculated uncertainties in 
the stiffness and damping coefficients for seals Si, HCl, HC2, and HC3 are given 
in Table 5. 


Table 5. Maximum Uncertainties. 


Seal 

k 

C 

K 

SI 

HCl 

HC2 

HC3 

31.2 N/mm 
6.5 N/mm 
9.0 N/mm 
9.4 N/mm 

0.044 N-s/mm 
0.036 N-s/mm 
0.050 N-s/mm 
0.038 N-s/mm 

20.4 N/mm 
16.9 N/mm 

23.4 N/mm 
18.0 N/mm 


Direct Damping 

All three seal analyses predict C well, especially for the honeycomb seals. For 
seal SI, predictions sire generally within 30% of the experimental results. For the 
honeycomb seals, most predictions are within 20% of the experimental results. 
Among the analyses, there is little difference in the predicted dependence of C on 
U 0 a and Pr. Plots of C versus u$ a and Pr are shown in figures 14-17. Plots of C 
versus rotor speed are shown in figures 18 and 19. For seal Si, there is no difference 
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Figure 18. C versus rotor speed for seals Si and HC1 of Table 1. 
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Figure 19. C versus rotor speed for seals HC2 and HC3 of Table 1 
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in the predicted rotor speed dependence of C for the three analyses. The fully- 
developed-flow and empirical-friction-factor analyses predict that C increases as 
w increases for all four test seals. For the honeycomb-stator seals, the flat-plate- 
flow analysis predicts that C decreases as the rotor speed is increased from 3000 
cpm to about 13000 cpm, and then increases as rotor speed is further increased 
to 16000 cpm. In general, the rotor speed dependence of C is best predicted by 
the fully-developed-flow and empirical-friction-factor analyses. 

Cross-Coupled Stiffness 

Typical plots of k versus rotor speed are shown in figures 20 and 21. The 
negative experimental croes-coupled stiffnesses in the figures may be due to sec- 
ondary flow patterns which are predicted by Tam et al. (1987) for smooth surfaces, 
low rotor speeds, and U$ a < 0, but are not modelled by bulk-flow analyses. For 
k > 0, the increase in k due to an increase in u is best predicted by the flat- 
plate-flow and empirical-friction-factor analyses. Due to higher rotor drag in the 
entrance region, the predictions of the flat-plate-flow and empirical-friction-factor 
analyses are more speed dependent (i.e., Ak/Au is greater) than those of the 
fully-developed-flow analysis. There is generally little difference in the predicted 
dependence of k on u for the flat-plate-flow and empirical-friction-factor analyses. 

Figures 22 and 23 show the dependence of k on u& a . The fully-developed-flow 
and flat-plate- flow analyses overpredict the increase in k due to am increase in u^ 0 , 
especially for the honeycomb seals. Because of the high /« near the seal inlet, the 
empirical-friction-factor analysis predicts the increase in k due to an increase in 
u ffa well. 

For seals Si, HCl, HC2, and HC3, figures 24 and 25 show k versus Pr for 
ue a « 1.0. For seals Si, HCl, and HC2, the empirical-friction-factor analysis 
predicts the dependence of k on Pr better than the fully-developed-flow and flat- 
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Figure 20. k versus rotor speed for seals SI and HCl of Table 1. 
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Figure 21. k versus rotor speed for seals HC2 and HC3 of Table 1 
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Figure 24. k versus P r for seals Si and HCI of Table 1. 
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plate-flow analyses. There is little difference between the predictions of the fully- 
developed-flow and flat-plate-flow analyses, and they predict the dependence of k 
on Pr for seal HC3 better than the empirical-friction-factor analysis. 

Direct Stiffness 

All three analyses predict K poorly, especially for the honeycomb seals. How- 
ever, the direct stiffness of a gas seal normally affects turbomachinery rotordy- 
namics less than k and C. Figures 26-31 show typical comparisons of predicted 
and experimental K for the test seals. A probable cause for the disagreement 
between experiment and analysis for the honeycomb seals is poor modeling of the 
honeycomb stator shear stresses. The friction factor characteristics of the honey- 
comb stators were obtained from tests of the seals at a clearance of 0.41 mm. The 
dependence of the friction characteristics of honeycomb surfaces on seal clearance 
has not been investigated. 
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Figure 27. K versus rotor speed for seals HC2 and HC3 of Table 1 
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Figure 28. K versus ue a for seals Si and HC1 of Table 1 



59 


SERL HC2, Pa = 3 .08 bars 
ROTOR 5PEED=9500 cpm 



.2 . B 1.0 


circumferential velocity rrtio 



Figure 29. K versus uo a for seals HC2 and HC3 of Table 1 
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Figure 31. K versus P r for seals HC2 and HC3 of Table 1 
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CHAPTER DC 
CONCLUSIONS 

Nelson’s (1984) rotordynamic analysis for centered, annular gas seals has 
been revised by including the fluid circumferential velocity in the Mach number 
definition. Two friction factor models have been developed and used in annular 
seal analyses. In the flat-plate-flow model, the friction factor in the seal entrance 
region is based on flat plate flow theory. The empirical-friction-factor model 
uses empirical entrance region and exit region friction factors. A method for 
determining the friction characteristics of honeycomb stators has been described. 
The results from static and dynamic tests of seals with smooth rotors and smooth 
or honeycomb stators have been compared to the predictions of: 

(a) a fully-developed-flow friction factor analysis; 

(b) the flat-plate-flow analysis; 

(c) the empirical-friction-factor analysis. 

The comparisons in the preceding chapters support the following conclusions: 

(a) Presently, only bulk-flow analyses can be used to predict the rotordynamic 
characteristics of honeycomb-stator seals. The three-dimensional flow analysis of 
Tam et al. (1987) provides a qualitative explanation for the inability of bulk-flow 
analyses to predict negative experimental k values for no fluid prerotation and 
low rotor speeds. However, the stator-wall boundary conditions assumed by Tam 
et al., and by Nordmann et al. (1987) in another three-dimensional seal analysis, 
make the use of their analyses questionable for modelling honeycomb-stator seals. 

(b) For a bulk-flow annular gas seal analysis, cross-coupled stiffness predic- 
tions are significantly improved by accurate modelling of the wall shear stresses 
in the seal. Compared to the fully-developed-flow and flat-plate-flow analyses, 
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the empirical-friction-factor analysis predicts k best for annular gas seals. Exper- 
imental rotor surface speeds are lower than those in many seal applications, so 
the improved ability of the flat-plate-flow and empirical-friction-factor analyses to 
predict rotor speed dependence is attractive. Since the inlet circumferential ve- 
locity ratio (Ufa) is between 0.5 and 1.0 in many seal applications, the improved 
ability of the empirical-friction-factor analysis to predict the dependence of k on 
u$ 0 is especially attractive. 

(c) Since the friction factor is high for a flow distance of about fifty 6eal 
clearances, the importance of using an entrance-region shear-stress model, instead 
of an abrupt entrance-loss, increases as the ratio of seal-clearance to seal-length 
{H/L) increases. 

(d) The fully-developed-flow and flat-plate-flow analyses overpredict the pres- 
sure loss in the last 10% to 25% of a seal. Generally, the fully-developed-flow 
and flat-plate-flow analyses underpredict the pressure loss at the seal inlet. The 
empirical-friction-factor analysis, with an exit region friction factor and a higher 
entrance region friction factor than the other two analyses, predicts seal pressure 
gradients well. The dependence on seal length and clearance of the empirical 
friction factor model developed in Chapter VI should be investigated further. 

(e) Further experimental investigation of the entrance-region shear stresses 
in seals would improve the present entrance-region friction model. The “forty- 
diameter rule” used in the flat-plate-flow and empirical-friction-factor analyses is 
based on the experimental results of Keenan and Neumann (1946) and Shapiro 
and Smith (1948), for which the minimum duct diameter was 9.5 mm. The results 
of Barbin and Jones (1963) and Wang and Tullis (1974) show that wall shear stress 
attains the fully developed flow value within fifteen diameters of the entrance of 
larger ducts (200 mm to 300 mm diameters). The experimental results of the 



64 


present study agree with those of Keenan and Neumann and Shapiro and Smith. 

(f) The honeycomb-stator friction characteristics derived by the method de- 
scribed in chapter V have been used effectively for predicting static and dynamic 
characteristics of honeycomb-stator /smooth-rotor seals. Only the direct stiffness 
is predicted poorly by the empirical-friction-factor analysis. However, the method 
described in Chapter V does not model the dependence of the friction charac- 
teristics of honeycomb surfaces on seal clearance. Flat-plate tests of honeycomb 
surfaces which include variation of the clearance between the plates will provide 
an improved honeycomb shear-stress model. Implementation of flat-plate test 
results in the present seal analysis should improve predictions of the dynamic 
characteristics of honeycomb-stator/smooth-rotor seals. 

(g) The entrance-region shear stress model used in the flat-plate-flow analysis 
is also applicable to incompressible seal analysis. The applicability to incompress- 
ible seal analysis of the friction model used in the empirical-friction-factor analysis 
must be determined experimentally. 

(h) All three analyses predict K poorly. For gas seals, this coefficient generally 
does not have a major influence on rotordynamics. 
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APPENDIX 

Uncertainty in Friction Factors 

The friction factor uncertainties were determined using the method described 
by Holman (1978). The uncertainty wr in a result R which is a function of n 
primary measurements Xi to x n with uncertainties u>i to xv n is 



The primary measurements in the friction factor calculations are seal clearance, 
pressure, flow rate, and stagnation temperature. The uncertainties in these mea- 
surements are about 13 ^m, 0.0055 bars, 1%, and 0.6°K, respectively. As a result, 
the estimated uncertainty in the Mach number calculation is about 1%. Friction 
is high near the seal inlet, where dM/dz is small. Using equation (46), a small 
uncertainty in dM/dz causes a large uncertainty in the friction factor calculation 
when / e is high and dM/dz is low. dM/dz is calculated using the curve fit of the 
experimental Mach numbers. For the smooth seals, the curve fit is within 1.5% of 
the experimental Mach numbers. For the honeycomb seals, the curve fit is within 
from 1% to 4% of the experimental Mach numbers. 

Accuracy of mo and no 

Figure A1 shows a comparison between average experimental friction factors 
(designated by the letter “x”) and the friction factors predicted (solid line) using 
equation (2) and the friction coefficients in table 4. The error bands show the 
range of experimental values used in calculating the average friction factors. 
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